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Abstract 

We consider the production of neutralino pairs x'iX'j ^ high energy hadron collider, 
putting a special emphasis on the case where one of them is the lightest neutralino Xiy 
possibly constituting the main Dark Matter component. At tree level, the only relevant 
subprocess is qq — > XiXj; while the subprocess gg 'XaX^ first appears at the one loop 
level. Explicit expressions for the gg-helicity amplitudes are presented, including the tree 
level contributions and the leading-log one loop radiative corrections. For the one-loop 
gg Xi'Xj process, a numerical code named PLATONggnn is released, allowing the 
computation of da/dt in any MSSM model with real soft breaking parameters. It turns 
out that acceptable MSSM benchmark models exist for which the qq and the gluonic 
contributions may give comparable effects at LHC, due to the enhanced gluonic structure 
functions at low fractional momenta. Depending on the values of the MSSM parameters, 
we find that the LHC neutralino pair production may provide sensitive tests of SUSY 
models generating neutralino Dark Matter. 

PACS numbers: 12.15.-y, 12.15.Lk, 13.75. Cs, 14.80.Ly 
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1 Introduction 



Neutralino production at hadron colliders is an important part of the program of Super- 
symmetry (SUSY) searches [H [2]. One special reason is related to the possibility that 
Xi, the lightest neutralino state, is in fact the Lightest Supersymmetric Particle (LSP) 
[3]. This has two particular consequences; the first concerning the supersymmetric spec- 
troscopy (chains and rates of decays) in R-parity conserving models [T]; while the second 
largely determining the search for Dark Matter (DM) |4]. 

DM detection in such a case is expected to occur either in a direct way {e.g. through 
the observation of nucleus recoil in Xi^ ~^ Xi^ elastic scattering); or in an indirect way, 
by observing modifications of the cosmic spectrum of particles like photons, positrons, 
antiprotons etc., due to contributions from XiXi annihilation [S]. Concerning the indi- 
rect way, we have presented in two previous papers the results of a complete one-loop 
computation for the processes XiXj ~^ 11)99 involved in DM annihilation [6j, as well 
as the results for the reversed process of neutralino pair production at a photon-photon 
collider [7]. In [H] we have also emphasized that in certain benchmark MSSM models 
the gluon-gluon channel may be important for determining the neutralino relic density 

We would expect therefore, that for neutralino pair production at a high energy hadron 
collider like LHC, kinematical domains may exist where the gluon structure function of 
the proton is so large pJ], that the one-loop gluon annihilation contribution may in 
fact be bigger than the tree level qq contribution. The precise study of such neutralino- 
pair production process at LHC, through the subprocesses qq Xi'Xj and gg XiXp 
constitutes the aim of the present paper. 

First, we present the hehcity amplitudes and cross sections of the subprocess qq 
Xi'Xj cit tree level. At this level, such a process has been studied long ago e.g. in p!2| [T3]. 
We go beyond this though by also exploring the Fermi statistics and CP (for real MSSM 
parameters) constraints on these amplitudes, which strongly reduce their number and also 
serve as a check of the calculations. 

In a second step, and in order to check the possible existence of important one-loop 
electroweak (EW) contributions to these qq amplitudes, the leading (a In^ s) and sub- 
leading {a In s) contributions are included, following the procedure established in [T^ \TE\ 
dn]. These EW corrections reduce the overall tree level magnitude of the amplitudes by an 
amount that can reach the few tens of percent level for the kinematical domain attainable 
at LHC. In the same direction acts also the SUSY QCD contributior0 calculated according 
to the rules established to order as In s in [161 [H] • We should emphasize at this point 
though, that these EW and SUSY QCD corrections should be considered in addition to 
the pure QCD (leading and next-to-leading) corrections which strongly increase the tree 
level amplitudes, as found in |T7]. 

We then turn to the one-loop subprocess gg X^x'j) ^oi which the helicity amplitudes 
are calculated using the set of diagrams established in [6l [7]. There, the neutralino 
annihilation amplitudes XiX'j 99 were calculated under any kinematical conditions; 

^ These involve QCD interactions explicitly affecting SUSY particle exchanges |16) . 
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but the accompanying numerical codes compute the neutrahno MSSM annihilation cross 
section to gluons only at the appropriate for dark matter threshold region [TS] . Using these 
results, the numerical code PLATONggnn has been also constructed, which calculates 
the reversed process cross section da{gg —>■ 'j<^X^)/dt for any (s,t)-values and any MSSM 
model with real soft breaking parameters |18j . 

We then compute the LHC cross sections for PP XiXj + •••5 by convoluting the 
gg and qq subprocess cross sections, with the corresponding quark and gluon distribution 
functions in the initial protons P. We then discuss the contributions of both subprocess to 
several observables (invariant mass, transverse momentum and angular distributions) and 
we give illustrations for an extensive set of benchmark models in MSSM. As we will see 
below, depending on the choice of MSSM parameters and the kinematical regions looked 
at, the one loop gg — > XiXj* subprocess may occasionally give comparable or even larger 
effects, than the tree level qq — ^ Xi'Xj' oiie. 

These results imply an interesting complementarity between the future LHC measure- 
ments, the related 77 Xi'Xj measurements at a future Linear Collider and the Dark 
Matter searches in cosmic experiments. 

The contents of the paper is the following. Sect. 2 is devoted to the process qq XiXj- 
The general properties of the helicity amplitudes are studied in the subsection 2.1, where 
the seven basic independent amplitudes are identified. The tree-level helicity amplitudes 
and cross sections are subsequently presented in Section 2.2 and Appendix A.l; while 
the electroweak and SUSY QCD corrections to the hehcity amplitudes, at leading and 
subleading logarithmic accuracy, are given in Section 2.3 and Appendix A. 2. In Sect. 3, the 
one loop process gg —>■ XiX^ is presented. Applications to neutralino pair production at 
LHC using the parton formalism are given in Sect. 4, where the numerical results are also 
discussed. The concluding remarks are given in Sect. 5, while the parton model kinematics 
is detailed in Appendix B. 

2 The subprocess qq XiX^j 

2.1 Generalities about Helicity amplitudes for qq X^X^ 

For an incoming gg-pair, and an outgoing pair of neutralinos, the process is written as 

Ai) q{q2, A2) X^iiPi, Ti) x°{Pj, Tj) , (1) 

where (qi, q2,Pi,Pj) ancl^ (Ai, A2, Tj, r,) are the momenta and hehcities of the incoming 
and outgoing particles. The usual Mandelstam variables for the subprocess are defined as 

s = (gi + 52)^ = {Pt+Pjf , i = {qi-pif = {q2-Pjf , u = {qi-pjf = {q2-Pif- (2) 

Since the top quark structure function is vanishing in the proton and the other quarks 
are not too heavy, the incoming q and g in ([1]) are taken as massless as far as the kinematics 

^The possible values of the hclicities Ai,A2, Ti,Tj are, as usually, taken as ±1/2. 
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is concerned, but we keep the potentially large (particularly for the third family) Yukawa 
contributions to the couplings. Finally (mj, rrij) denote the (x^, Xj) masses, respectively. 
The hehcity amplitude for process ([T]) is denoted as 

^Zxr,nr,m , (3) 

where 6* is scattering angle in the cm. of the subprocess. In these amplitudes, q and Xj 
are treated as particles No2 in the JW conventions ^19j . To consistently take into account 
the Majorana nature of the neutralinos, we always describe the Nol neutralino x° through 
a positive energy Dirac wave function, while the No2 JW particle Xj is described through 
a negative energy on^. Fermi statistics for the final neutralinos then implied 

^^a,.^(^*) = (-1)'^"'^^:a.^.(^-^*) ' (4) 
while CP-invariance, valid for real soft breaking and /i parameters, gives 

where rji = ±1 is the CP-eigenvalue of the Xi-neutralino [2T] . 

On the basis of (jU [5]) all gg-amplitudes may be expressed in terms of seven basic ones , 
selected as 



The other amplitudes are determined from these through 
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We also note that (jU [5]) imply the relations 

■^The same convention is also followed for the gg XiX^ treated below. 

^We note that (|3|), which is induced by the anticommuting nature of the Fermionic fields, does not 
generally agree with the neutralino (anti) symmetry property assumed in [20] . 
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In terms of these helicity amplitudes, the unpolarized differential subprocess cross 
section may expressed as 



di 1927rsn2 
where the kinematics are defined in ([2]) and flB.llHB.TGj) . 



\^Xi,\2;n,Tj\^ , (9) 



2.2 Born amplitudes and cross sections. 

The Born amplitude for the process in ([T]) contains three diagrams (see FigHKbc) involving 
s, t and u channel exchanges and written as [13] 

ipVB _ pijBs _|_ pijBt _|_ pijBu qn 



mZ.r, = - .2,2^s .A ^h'[9,LPL + g.RPRHq) ■ Ua,[N'^ Pl ' N^^*PRh, 

^ ^-^ t — mi 

n " 

nZ.r, = - + ^^*i^n)PLK ■ v-{Af{qn)PL + Af{q^)PnHq), (11) 

n " 

where the index n refers to the summation over the exchanged L- and R- squarks of the 
same flavor in the t- and u-channel, Pl/r = (1 T 75)/2, and {i, j) describe the final 
neutralinos. 

Explicit expressions for the seven basic helicity amplitudes listed in (|6]), are given in 
(lAH ICT 1X31 \EM in Appendix A.l. 

• They involve the L and R Zgg-couplings defined as 

Zqq -^^^;^'^3,LPL + 9gRPR] , (12) 

where 

g,L = 2Pgil-2sl,\Q,\) , 9,R=-2sl,Q, , (13) 
with Ig, Qq being the isospin and charge of the various g^-quarks. 
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• The Z-neutralino couplings satisfying 

ZX'X' ^ TT^^'h.LPL + gi.BPB] (14) 
ISwCw 

with 

g,,L = N'^ ^ ZtZ^j - Z^:Z^^ , g,,n = -N''* , (15) 
where denotes the neutrahno mixing matrix in the notation of [22]). 

• And the neutrahno-quark-squark couphngs 

qqnt Af{qn)PL + A^{qn)PR, (16) 

where 

^fi^n) = , A^iu,) = - . Z^ 

^ V2Mvi/Svp'Sm/3 

= -^Z-C , Af{d,) = -Z^,: (IT) 

^ v2Mvi/Stycos/^ 

In (|T7l) . (g = M, d) refer to the incoming up and down quark (antiquark) of any 
familjcl, while (g„ = gz., q^) denote the corresponding squarks. We also note that the 
mixing matrices Z^ in ( TT5l [T71) . control the Bino, Wino, Higgsino components of the 
neutralino in the Zx^x^ and qqx^ coupling [22l [21]. Finally, we remark that s-channel 
Born part F^^^^^_^^ gives non- vanishing contributions only for purely higgsino production, 
whereas the t, w-channel Born parts for purely gaugino. 

One can then compute the differential cross section either through ^ and the helicity 
amplitudes in flA.H IA.2[ IA.31 IA.4p . or directly by the trace procedure giving 



da{qq ^ Xlx]) _ 1 ^1^'^ , 

' 'tu\ 



, „r. -9(0) V-ss + ht + luu ~ '^Ist + 2/su — 2/j 
dt 1927rs'' V2/ 

wher^ 

Iss = , 'y})^ 'Hj ^'^N-ni^l - i){m] - i) + (m^ - u){m^ - u)] 



^As usual, we will only consider non- vanishing structure functions for incoming u, d, s, c and b quarks. 
^Analogous expression for gaugino production in e^e+ Colliders have appeared in [23] . 
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+Af{q,)Af*m , 

{m] -u){m] -u)^ ^^^^^ ^ ^^^^ 



1st 
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+Af*iq,)Afiqi)] , 



[Ar{qk)Af{qi) + At{q,)A^{mA^*i^k)Af{q, 



Re\ [N^^*g,nA^iq,)Af*{q,) - N^^g,LAi^iq,)Ai*iq,)] 
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'm] - i){m] - i) + [iV*^*(7,z.Af - N^' g,RAf{q^Af{q^)\m,mis 



{m\ - u){m] -u) + [N'^*g,nAf*{qk)Af{qk) - N'^ g.LAf* {qk)Af{qk)]m,m,s 



Ltu 



{u-mt){t-ml) 



Re\[A^{qu)A^{qi)A'^\q,)A'^\qi) 



+Af{q,)Af{qi)Af\qu)Af\qi)]m,mfs + -^^^^ 

+Af*{q,)A^*{qi)A^{q,)Af{q,)\ [K - " i) + K - " ^) 

-s{s-m^^-m])]^ , (19) 

where m^, mj are the neutrahno masses and A^*-^ are defined in (fT5i) . 

The resuhs in eq.(|T8l) disagree with those of [20], where symmetry properties of the 
XiXj states which are different from those in (jl]) have been used. 



2.3 One loop electroweak and SUSY QCD corrections 

to qq 

In principle, one loop EW corrections for qq — > XiX? should be taken into account, partic- 
ularly because the energy reach at LHC is so big, that the large logarithmic contributions 
to the amplitudes may reach the few tens of percent level [HI |T5l [16] . In the models we 
have considered, this implies a reduction of the cross sections sometimes by almost a fac- 
tor of two, while preserving their shape. Since the non-logarithmic one-loop contributions 
seem to lie at the few percent level, which is also the level of the expected experimental 
accuracy, it may be adequate to ignore these difficult to calculate effects in {qq — > Xi'Xj) 
at LHC energies. 

In this section we present therefore, the leading and subleading EW logarithmic con- 
tributions to the qq helicity amplitudes, following [HI [151 [16], where applications 
for LC and LHC have been given. They are separated into three types of terms which 
are: 
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Universal electroweak (EW) terms. These are process-independent terms appearing 
as correction factors to the Born amphtude. They consist of "gauge" and "Yukawa" 
contributions associated to each external line and determined by its quantum num- 
bers and chirality. Their expressions for a quark or neutralino line are respectively 
determined as follows: 

External quark line of chirality a = R : Since all quarks are taken as massless 
as far as the kinematics are concernecill, the quark lines correspond to a definite 
chirality a. The induced correction then is 



(20) 



where F^^^ describes the corresponding Born amplitude from (ITOl) involving an 
external quark line of chirality a (which at high energies is essentially equivalent 
to the helicity), while (i,j)-count the mass eigenstates of the neutralinos. The 
coefficient in f l20p is written as 



c{q)a = c{q, gauge)a + c{q, yuk), 
with the gauge contribution being 

r/„(/o + 



c(g, gauge)a 



a 
8^ 
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w 



m 



w 



(21) 



(22) 



where Ig is the full isospin of the quark with chirality a = L or R, and Yg = 

'^{Qq ~ H ) defines its hypercharge in terms of the third isospin component Iq . 
Correspondingly, the Yukawa term (for b, t quarks only) is 



c(g, yuk)a 



a 



In- 
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w 



+ ■ 



rriT 
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mt 



sin^ (3 rri^^ cos^ [3. 

2 



.m%^ cos^ (3 ^'i^' ' ~'~ mir sin^ [3 



4 smM 4'\V2, 



5aR 



(23) 



We note that an external antiquark line should be counted separately giving an ad- 
ditional contribution determined by the same formulae fl20ti23p . Moreover, the same 
formulae describe also the logarithmic contributions associated with each external 
squark, anti-squark, lepton or slepton line [T H \Wl \W\ . 

For an external neutralino line of chirality b, it is convenient to use a matrix 



notation 



(24) 



and to separate the higgsino from the gaugino components of the matrix elements: 



c{Xi Xi')b = c{xi Xv higgsino, gauge)?, + c{xi Xu higgsino, yuk)^ 
+ c(xl'x° gaugino, gauge)6 , 



(25) 



^The large third family Yukawa terms appear only as couplings and do not concern the kinematics. 
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with 



c(x?X° higgsino, gauge)b 
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ciXiXi' higgsino, yuk)b = - 
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(26) 



The logarithmic contributions associated to the purely higgsino s-channel Born am- 
plitude F^''^^ o^ FigH^ only involve the "higgsino, gauge" and "higgsino, Yukawa" 
elements, whereas the contributions associated to the purely gaugino t- and u- Born 
amplitudes i?*^^*'""^ only involve the "gaugino, gauge" elements. 

Angular and process dependent terms. They originate from diagrams involving W 
internal lines supplying soft-infrared In^ t or In^ u terms. Diagrams with internal Z 
lines are negligible, since their contributions turn out to be orthogonal to the Born 
terms and cannot interfere with them. The contributing diagrams therefore consist 
of boxes with an intermediate WW pair in the s-channel, triangles involving a single 
W connected to a squark exchange in the t or u channels, and boxes involving a 
single W and a squark in the t or u channels. 

Renormalization Group (RG) terms. They arise from intermediate Z boson Born 
terms contributing to Higgsino production only and inducing running effects to the 
corresponding g and g' gauge couplings; see Fig. la. In terms of the s-channel Born 
amplitudes of ffTOl) . they are written as 



pRG 



1 



47r2 



where g^ = e^/s^, g'^ = e^/c 



~ 3 Ug Uh 



dg^ 
1 

~ 4 



~.dF'^^' 



dg 



12 



5 

6 ^ 



11 

T 



(27) 



(2^ 



and Ca = "2, Ug = 3, = 2 in MSSM. Applying this procedure to the 7 basic Born 
helicity amplitudes of Appendix A.l, implies the substitutions 



„2 2 

e^gqR 
„2 2 



47r2 

(2gg 



In- 



/4| 



i-2|g,|]} 



47r2 



iPW'} 



(29) 



^See also RUh. 
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The resulting expression for the seven basic hehcity amphtudes of (jS]), are given in 
Appendix A. 2. For what concerns the magnitude of the various corrections, the following 
general comments can be made. In the LHC domain (say for ~ 1 TeV) a single 
ln(s/m^) and a squared ln^(s/m^) give enhancement factors of about 5 and 25, re- 
spectively. So one expects that the corrections are of the order of —5% for the quark 
or higgsino-gauge terms, —25% for the gaugino gauge ones, and —10% for the higgsino- 
Yukawa terms (depending on tan/3 value). The angular dependent terms have a more 
complicated structure; their sign is not fixed, while their magnitude can reach the 10% 
level. The addition of these various electroweak terms is strongly model dependent, espe- 
cially due to the Zj^ matrix elements controlling the amount of the higgsino and gaugino 
components of the neutralinos. The net EW effect on the amplitude though, is essentially 
always negative and can easily reach the several tens of percent level. 

The logarithmic SUSY QCD corrections for quark-antiquark processes at order as are 
given by [17, 16j, 

(30) 



pSUSY QCD ^ pBorn 

Xq,Xq;Ti,T j Xq.Xq', 



This simple logarithmic terms arise from diagrams containing virtual squarks and gluinos 
interacting via SUSY QCD couplings; such diagrams produce no In^ s-terms. By itself this 
single log term is of course negative, and it would remain around —5% in the observable 
LHC domain. In addition to these though, we should always also consider the pure 
QCD leading and next-to- leading mass dependent corrections which, as shown in ref. |17j. 
result into a K factor that turns out to be positive and of the order of 30%. These later 
corrections are not discussed in the present paper. 



3 The one loop process gg XiXj 

This process first appears at the one loop level through the triangle and box diagrams 
fully listed in |6]. These diagrams basically involve gluon-quark-squark and neutralino- 
chargino-squark couplings; no gluino can appear at this order. Accidental degeneracies 
between the neutralino masses and squark masses can give some enhancement effects. In 
addition, single Z, h^, or exchanges in the s-channel, can also give enhancements 
and resonance effects at the corresponding cm. energies. These situations are rather 
similar to those already mentioned for the 77 — ^ XiXj process in [7]. 
The helicity amplitudes for the process 

giQl, /^l)^(92, /i2) X^iiPi, ^»)X?(Pj, Tj) , (31) 

are denoted as 

^...^(^*) ^ (32) 

where the momenta and helicities of the incoming gluons and outgoing neutralinos are 
defined, and 9* again denotes the cm. scattering angle. We use the same (x^, X?) 
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conventions as in the qq case and in implying 

^,,.;^,^(^*) = (-1)^-^^ - (33) 

from Xi'Xj fermion-antisymmetry, and 

from (7(7-boson symmetry. 

If the MSSM breaking parameters and the Higgs parameter /i are real, then CP in- 
variance holds, implying 

F%,,-,,,-r.,-r,m = (-l)--^-(---^) r^.r^,i^,,,.^,^(r) , (35) 

where 77^, 77^ = ±1 are the CP-eigenvalues of the two produced neutralino^. In such a 



case, time inversion invariance implies the same helicity amplitudes for the process (1311) 
and its inverse. Combining (l33l IMl l35ll . we get 

Fii,^,nr,m = (-i)^^-^^^^-"i^^,,;^ = ^.^.n„_,,_^,_,(r) , (36) 

where the first part comes from ( 133| IMl) alone, while for the last part the CP-invariance 
relation (l35l) is also used. 

In terms of these helicity amplitudes, the unpolarized differential subprocess cross 
section is 



(2) 5Z I^Mi/^2r,r,r- (37) 



di 40967rs2 

Together with the present paper, we release in [18] the numerical code PLATONggnn, 
which calculates the differential cross section (1371) as a functions of 9* and s, for any set 
of real SUSY parameters at the electroweak scale. 



The above amplitudes are basically of order as/vr weaker than the tree level qq am- 
plitudes of the preceding Section. In certain SUSY models though, this reduction can be 
partially compensated by the aforementioned enhancement factors. But in practice, the 
most important feature at LHC is the relative size of gg and qq distribution functions 
inside the proton; i.e. the fact that at "low" subenergies the gg fluxes are much larger 
than qq ones. Because of this, and as we see in Section 4, there are benchmark models 
where the gg contribution to neutralino-neutralino production at LHC, is larger than the 
qq one. 



4 XiXj distributions in Proton-Proton collisions 

In this section we discuss numerical results for the process PP X^Xj + .... at LHC (cm. 
energy ^/s = 14 TeV) generated by the subprocesses gg, qq 'XaX^i {'Xa is always taken 

^We follow the same notation as in e.g. l21] . 
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heavier than All necessary formulae describing the kinematics for two massive final 
particles are presented in Appendix B. These allow the computation of the neutralino- 
pair invariant mass distributions da/ds, the transverse energy distribution of the heavier 
neutralino da/dxri, and the angular distribution in the neutralino-neutralino center of 



mass da/dxi described through Xi defined in (16.19^ : compare the formalisms in Sections 
B.2.3, B.2.1 and B.2.4 respectively. 

The present formalism allows of course to compute any x^x^ channel for any MSSM 
model with real soft breaking and fi parameters, using code PLATONggnn released in P^. 
In this paper we also assumed that the Xi escapes the detector without being observed^"!, 
so that the identification of e.g. the x^Xi-pi'oduction is only done through the detection 
of X2- Our illustrations are restricted to the X2X1 X2X2 channels. 

As an example, for the quark and gluon distribution functions inside the proton we 
use the MRST2003c package [TTj at the scale (lB.30p . Using this and the 31 benchmarks 
models [8],[9l|10] already considered in the previous paper^^ [6l [7] , we have made numerical 
computations of the above three single variable distributions. Two main features come 
out from this study. 

At high invariant subprocess energies (close to 1 TeV), the gg contribution becomes 
negligible compared to the gg-ones. This is due to two effects; the gg subprocess cross 
section is reduced by the one loop factor Os/tt compared to the tree level gg-subprocess, 
while the gg luminosity is comparable to (or even weaker than) the gg one. 

On the opposite energy site, within a few hundreds of GeV above threshold, the gg 
flux may be so large, that the gg contribution may compete or even overpass the gg- 
contribution by a factor of 10 or more in some of the benchmark models. This is further 
enhanced in cases where the or Higgs boson can couple to X2X1 o^' X2X2 channels. 

For the illustrations presented here, we have mainly selected those of the benchmark 
models mentioned above where the competition between the gg and gg contributions is 
most spectacular. Thus, in Figsl2ll5] we show the invariant mass distribution da/ds in 
the cases where there is a resonant enhancement in the gg contribution (models SPSlal, 
SPSS, SPS6, SPSS ,[8J), while in FigsEl El the 99 contribution is smooth but irnportant 
at low subenergies (models SPS7, CDG24, [U [ID]), and it is further reduced iio SPS4 

(FiglHD. 

In all cases, one sees that the gg contribution to da/ds has a larger slope than the gg 
contribution; the effect being mainly due to the behavior of the gluon distribution func- 
tions. The precise magnitude of the gg contribution is however strongly model dependent 
and arises as the result of many features of the SUSY spectrum involved in the one loop 
diagrams contributing to gg XiX'j [S]- 

As already discussed at the end of Sect. 2. 3, the one loop logarithmic electroweak 
corrections to gg XiX'j tend to reduce the size of the gg tree level cross sections; see 
also [m [151 [iGj . They are strongly model dependent (due to the neutralino mixing matrix 



-'^^This would be the case if Xi is the stable LSP. 

^^The high scale values of the defining parameters of these models are listed in Tables 1,2,3 of [6]. 
^^The gg contribution may become even smaller in some of the other benchmark models mentioned 
above. 
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and tan f3), and can reach several tens percent for the amphtudes. For certain models, the 
addition of the various terms can lead to a reduction of the size of the qq cross sections by 
almost a factor two, but they do not strongly modify their shapes. This reduction will of 
course be somewhat reduced by the enhancements induced by the pure QCD effects pTTj . 
but it is nevertheless sufficient to increase the relative importance of the gg contribution. 

The features observed in the da/ds distributions of Figs|2ll8l can also be seen in the 
transverse energy distributions da/dxTi- In order to not multiply the number of figures, 
we only give illustrations of this fact for two typical models, SPSla in Fig.9 (where there 
is a resonance), and SPS7 in Fig. 10 (where there is no resonance). 

Finally we have examined the distribution da/dxi, which is essentially controlled by 
the neutralino-neutralino center of mass angular distribution; i. e. the cos 6* dependence 
discussed in Appendix B2.4 and ( IB. 191) . Such Xi distributions can also be used as a com- 
plementary test of the dynamics responsible for neutralino-neutralino production. Typical 
illustrations for models SPSla and SPS7 appear in Fig. 11 and Fig.l2 respectively. Con- 
cerning them we should remark, that the EW correction to the qq subprocess are reliable 
mainly in large Xi region, where the gg subprocess may also be important, especially if 
there is an or resonance effect. This appears as a threshold effect corresponding 
to the value of Xi above which the mass of the resonance lies within the allowed integra- 
tion domain for (]B.43|) . On the contrary, the large correction in the small Xi range of 
Figs. 11,12 should not be taken too seriously, since it is caused from a region where \t\ is 
small and the leading-log predictions not valid. 



5 Final discussion 

In this paper we have considered the neutralino pair production processes in proton-proton 
collisions at LHC. In the description we have taken into account the subprocess qq XaX^j 
calculated at the Born level as a first option, and as a second option we have included 
also leading and subleading logarithmic corrections. The genuine one loop gg — > 
subprocess, is fully taken into account. The description applies to any MSSM model with 
real soft breaking and fi parameters. Analytic expressions for the helicity amplitudes have 
been explicitly written for gg-subprocess, while a numerical code PLATONggnn is released 
allowing the computation of the rather involved da{gg — > XiX^)/dt, for any neutralino 
pair [IS]. 

After convoluting the subprocess cross sections with parton distribution functions, 
several observable distributions in PP XiX'j + ■■■■ at LHC, have been studied. For the 
applications, we have restricted to X2X1 and X2X2 production in the context of 31 bench- 
marks MSSM models also considered in [HlEj. A strong model dependence is observed, to 
which almost all aspects of the MSSM spectrum contribute through masses and mixing 
matrix elements. 

One of the most striking feature we have found is the important role of the gg sub- 
processes, which, although basically suppressed by the one loop as/71 factor, it may oc- 
casionally supply a larger contribution than the qq subprocess. This may occur close 
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to and slightly above threshold, where the gg luminosity could be sufficiently large to 
compensate for the Os/vr factor. In some models, the one loop gg amplitudes may be 
further enhanced by the presence of or Higgs boson resonances, and possibly also 
by accidental degeneracies between the neutralino and squark masses. 

We have also given detail illustrations for the invariant mass, transverse energy and 
angular distributions, for the case of seven benchmark models where the gg contribution is 
generally spectacular 0121110]; i-e. appearance of peaks, threshold effects etc. In the few 
hundred GeV subenergy domain, such structures of the da/ds distributions, appear in the 
range of 1 to 100 fb/TeV~^ and should be observable at LHC This may also be true for 
the da/dxTi and da/dxi distributions; compare Figs. 9-11. We should remember though, 
that there exist benchmark models also, where the gg contribution is rather marginal, as 
e.g. in model SPS4 [8j and others P [10]. 

In all cases, the effect of the one loop logarithmic corrections to the qq cross sections 
appear to be at the few tens of percent level or more, compared to the tree contribution, 
and should be taken into account in LHC computations. 

These features make the neutralino pair production processes rather interesting for 
testing the SUSY dynamics at LHC. The reason is that they provide tests which will 
be complementary to those addressing the cascade decays of initially produced colored 
SUSY particles to eventually Xi; which is here assumed to be the LSP; e.g. studies 
of mass spectra and decay branching ratios [2]. In particular, consistency checks should 
thus become available, allowing the strengthening of possible constraints on the validity 
of specific models. 

Moreover, in such neutralino pair production, the role of the Majorana nature of the 
final state particles is more prominent than in decays involving just one neutralino at a 
time. Since no such states (except possible the neutrinos) have been observed in the past, 
it would be interesting to have eventually some experimental support of our understanding 
of the Majorana nature. 

If Xi turns out to be an important or dominant component Dark Matter, the present 
calculation^, should also help in providing LHC constraints on the direct or indirect 
observations of the Dark Matter properties, whenever they will become available. 



For the same reason, Linear Collider studies should also be helpful [7]. 
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Appendix A 



A.l Tree level helicity amplitudes for qq xlx 



Using the notation of ( fTOl) . the Born contributions arising from the s-channel diagram in 
FigH^, to the seven basic hehcity amphtudes hsted in ([6]) consist of 



rpijBs 



pijBs 



pijBs 



-lijBs 



e^QqLVS 



- ml) 



(1 + cosr) Js - {rui - rUj^iN'^ - N'^*) 



e gqRVs 



(1 + cosr) Js - {rrii - rrij^iN'^ - N'^*) 



-sine* 



+ \/s - {mi + mjf{mi - mj){N'^ + A^*^*) 

rrvSin 9 



e^gqR 



+ Ws - {mi + mjf{mi - mj){N'3 + iV*^*) 



while 



■^ijBs 







(A.l) 



(A.2) 



for all (xj, rj)-valueo. Here 9* is the scattering angle in the cm. of the subprocess, and 
denotes the neutralino mixing matrix in the notation of [22] . 

The Born contributions to the seven basic helicity amplitudes of (E]), arising from the 
t-channel diagram in FigJTjD, are 



?ijBt 



TpijBt 



-^[\/ s - {rrii - rajy - ^Jrs-{m^ + m^f 

■ E iT^^^i^n)A^*i^n)] (1 + COSr) , 

„ yt m^ J 

- \ls-{m, + mjY 

■ (^—A^{qr.)Af*{q^)) (1 + cosr) , 



*In fact, the s-channel Born contributions vanish for all amplitudes with equal incoming helicities. 
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rpijBt 



rpijBt 



— -|^(mj + rrij)^ s — (rrij — rrijY — {rrii — rrij) ^J s — {rrii + rrij)'^ 



-|^(mj + rrij)^ s - [rrii - "^iJ + - "^Ja/s 



^ 1^ y s - (mi - mj)2 - - {mi + nijY 



pijBt 



pijBt 



- (mj + mj) W s - (mi - mj)^ + (m^ - m^) \/ s - (m^ + rrij) 



•Efr^A'^a)Af(gn)) sinr 



(A.3) 



while the corresponding contributions by the u-channel diagram in FigHb, are 



pijBu 



rpijBu 



• E (,-^Af{q^)Af*{q^)) (1 + cos^) 

-^l^y^s - {rrii -rrijY + ^Js - {rrii + rnj)^ 
■ E (^-^Afi^n)Af*{q„)) (1 + cosr) , 

- |^(mj + mj)\J s — {rrii — m^y + {rrii — ^j) \J ^ ~ + 

.J2(-l^^Af{qr.)Af*{q^)) sinr , 
\u — mi I 

— -|^(mj + mj)^ s — {mi — mjY — {m^ — mj) \J s — {mi + mj 

■E f-^^f(^~")^f*(^~'^)) ^i^^* ' 
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TpijBu 



s — {rrii — rrijY ~ y ^ ~ + ^i)^ 

\u — mi I 
^J^s - (mi -mjY + y^s - {rrii + rrijy 

^-^ \u — mi ■' I 

- \{mi + - {irii - mjY + {mi - mj) \/ s - {mi + mj^"^ 



E (-^Af{q.Mf*{q^)) sine* 



(A.4) 



The summation in flA.3t IA.4I) runs over the left and right squarks with the same flavor as 
the incoming quarks and mass m„. The couplings in (IA.lt IA.31 IA.4p have been defined in 



A. 2 Leading log helicity amplitudes for qq ')qx 



0.-0 

J' 



In this subsection we include the one-loop leading log contributions to the seven basic 
amplitudes of ([6]). We use the various couplings defined in (11 2111 7p and 



cos /3 



sin /? 



7N 



ryN — ryN 

5 _3 rig rih _ 1 



e/sw, g' = e/cw 

Pn = 

6^8 
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T 



while an In-symbol standing alone should be understood as 



In — > In- 



m 



(A.5) 
(A.6) 
(A.7) 

(A.8) 



w 



Denoting then by {q = u, d) the quark occurring in the initial state, and by (g' = (i, u) 
the corresponding companion quark belonging to the same SU(2) doublet, and describing 
by Ms the effective average mass for the the squarks of the same flavor as the incoming 
quarks, the seven basic amplitudes of ([6]) are written as 



a [In] 



m. 



'ail 



1447rs^c^ 



26c'v^)(21n-ln2) 



— In- 

37r 



m 



w 



■ + 



mt 



{5gt + 5gb) 



-+- 



ail + 
- levrsLcS 



Wj 



(21n-ln2) 



17 



3a' 



AO. 9 / 



+ COS 9" 



m|) 



2 



2 

(ryN*ryN ryNryN*\ ^^t 



- (mi 



2 

i + rrij) \Zj^^ Zj^j + ^4i^4j I - 2 o ~ V^3i + ^3i ■ 

olil ID 

_^^(3)a^[ln^] V^l (1 + cos^^^ ' 



2 

7iV yN*\_['];b 



t 



)srU / ^ 



(l + cosr) p r — ^ /T . ^ .2] 

/f + cos r ) I [Js-(m,-m,-)2 + Js-{mi + mjy 



) + ^2^(^i?^H^ + 6/fZ^ 



.[i2Zi 



'2i ^2j + 



+ 



{rrii 



■({2z^,zZ* + z^,z^;) _ 



- {rrii + mj)2 
\^^]-{2Z^,ZZ* 



1 



s - 



In- 



s - {rrii - TTij^ - - {rrii + 

{Y^^i^2i{zS*sw + ^if^z^;cw) + z^;{z^,sw + e/f )z2Tc^)] 



2 — 1 



+ 



I - 

s . 



4 2 4 ^q'i^q'j 

(1 + COS/?* 



'2^^^) + Z^i{Z^*sw - 6/^)^2^^^)] 



- rrijy + Js- {mi + m^y] 



■sw + Qlf^Z^*cw) + Z^C{Z^^sw + Qlf^Z^^cw)] 
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a 



s J \12s^cw 



4 ^q'j^q'i 



+ -^(i-2|g, 



J (s - m|) 



(1 + cosr)[(A^*^ - N'^*)\ s - {mi - m 



+{N'' + N'^*) Js - {mi + mjY)] [In] , 



(A.9) 



Q;[ln] 



++^ 
2 



■^(l±l^(21n-ln^)- ^In^ 
. UAns^c^ ' Stt mlr 



mt 



(21n-ln2)] 



tN i7N*\ 



[In] { (m, + m,) ^5 - (m, - m,)^ [(Z.fZ,^ - Z^, Z^.*)- , ^ 



cos^ pJ 



mi 



7N r7N*\ 



olli l) 



-{Z,, Z,^ + Z,,Z,^ )—^^ 



]} 



^^^^^ 12stcl -^^{ql} '^^' ^ ^^"^ - {mi-mjf - {mi - mj) ^ s - {mi + mj^] 

■[Z^nZ^rSw + 6/f )Z,^c;^) + Z^^{ZZ*s^ + 67f )zj*c^)] 
sin 0* 



u-m?{q )^^"^' mj)y^s - (m^ - mj)^ + (m^ - mj)^Js- {mi + mj)^] 
•[^2?(^iV + 3(2/f ))Z,^c^) + ^^(^if «^ + 3(27f c^)]} 
—if^^r-i — sin 6* < [(mj + mj) y s — {mi — m^Y + (m^ — m^) J s — {mi + m^)^] 

+[(mi + mj)^s - {mi - m^Y - {mi - mj)^s - {mi + m^Y] 



{2Z^X7 + zlz:*) 



)} 



s . 



{mi + mjY 



In — 
s . 



3m^cw 



m 



w 
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+ 



s J 

sin 6** 

u- Ml 



m 



W 



[{mi + TJij) \J s - {rrii — rrijY + (m^ - rrij)^ s — {rrii + rrijY] 



In- 



-u 



+ 



s 

S . 



Z^j{Z^*sw + Qlf^Z^*cw) + Z^*{Z^.sw + QlfZ^.cw) + 



m 
3ml, cw 



1 ^qi "^qj 

w 



] 



% + ^(l-2|Q,|) 



[s — mi 



+ {N'' + N''*){mi - mj)Js - (m^ + m,)2] [In] , 



(A.IO) 



(21n-ln^) 



«[ln] 



-Oat H 



47rs^m^ Vsin^ /3 ^ cos^/? 



etc , s 
^In 



(2 In -In') 



Sa^gqR sin ^* 



[In] (m, + mj)Js - (m, - m,)^ (Z.^Z,^ - Z^^Z^^*)^ 



sin'/? 



ryN 

■3j 



7N ryN*\ ^6 



COS^ 



^.r.a^mqmqi ^^^^^ I — 5 r/ , 



sinr[ln]f In^ 



In^ / 

^ — [{mi + mj)J s - {mi - mj^ 



u-Ml 



[m. 



- m,)^s-{mi + mjmZ^iZ^; + Z^-Zf )} 



+ ( AT^i + AT^i* ) {rui - mj ) ^s- (771^ + 771^)2] [i^] ^ 



(A.ll) 



-^In- 



.{on 
+ F|^;^+l(r) 



.47rc' 



(2 In -In') - 



a [In] 



w 



m: 



mt 



sin'/''*^^^^^ 



+ 2c^) 



167rs^c^ 



(2 In -In' 
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3Q;^fl'5i?\/^(l + cos 9*) 
16s^c^m^(s - m|) 



[In]^ - (m. - [{Z^Z^, - ^4^^^)-^ 



m 



.2 



( ryN* rvN rvNrvN*^ 

2 

-[Z^i ^3i + ^3^^3,- )^^^ 



(mi 



sin (5 



mt 



_(7V^i* + N'^)^s - {rrii + mj)2)][ln] , (A. 12) 



a [In] 



c^(l + 26c"^; 
2887rs2^c2^ 



+ 



a 



727rc 



w 



-^ln4 



37r" 



Stts 



sin^ (3 cos^ 



zi*z^^*{i-cose*) 



'(1 + cosr) 



Z2^*z,t(i + cosr) 



i - fh?{qL) 



t-Ml 



z^*{zZ*sw-Qif^z^;cw)\n-^ 



+ {6cw{z,fz^: + z^X) + ^Mz^*z^* + z^X)) 
(1 - COS e* 



u-Ml 



'zz\z^:s^-Qif^z^:c^)\n4+ 



qj ' "^ij qi 
-t 

s 

■N* ryN* , ryN ryN 



[^cw{Z^:Z^; + Z^X) + ^^w{ZZ*Zf + Z^X)) In- 



-u 



(A.13) 



_{e*) = 



a(l + 26c2 



2887rs2j,c2^, 



+ 



Lo7^c^r ' 



727rc2^ 



-^In. 



a 



[In] 



37r 



2v^[ln2] 



mi 



87is'^m,ly Vsin^ /? cos^ 



'Jqb 



a 



[In] 



sin^ /9 cos^ P 



{J s-(mi- mj)2 + Js-(mi + mj)'^] 



-z,yzg-(i-cosr) 

i-m?{qL) 
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_/(3) 



(1 -cosr 



t-Ml 



+ [6cw{Z^;Z^^ + Z^,Z^\) + Asw{Z^:Z^* + In- 



(1 + cosr) 



- Qlf^Z^Ccw) In— + 



6cH^(^2t^,t + ^2^^^) + ^s^{ZZ*Z^: + ZlIZ,^) ) In- 



-u 



(A. 14) 
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+ 



tts , s 
m 



:dot H ^„ 



8ad + 



V (21n-ln2 



Svr Sus^m^ Vsin^ /5 ' cos^ /? 



+ 



a [In] 



mr 



+ 



IGvrs^m^ Vsin^ /9 



+ 5qb) 



+ (mj - ■mj)J s - {rrii + ""^i)^] 



^mg[ln] sin 6* 



./(3) 



+ {mi - rrij) \ s - {rrii + "^i)^] 



[{rrii + ^j)\Js - {rrii - 



t-Ml 



Z^:{Z^*sw-Qlf^Z^*cw)\^^ 



^qi K^lj 



+ (6cw{z^;z^: + z^^z^^) + As^{z^:z^; + Zp^)) In^' 



u-Ml 



-t 



zS*{z^;s^-6i('^z^;c^)\n^ + 



{Qc^iZtZ^; + Z^X) + ^sw{Z^;Z^: + Z^,Z^^)) In-^] } 



(A. 15) 



expressed in terms of the Born amplitudes -^AiA2tt s-channel part F^^j^'^r r ^P" 

pearing 6.5-. in ffTOl [OlA.4p . 
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Appendix B 



B.l Parton model kinematics for x^'Xj production. 



The basic parton model expression for the hadron-hadron coUision A{qi)B{q2) — > X°(Pi) + 
^/ /^^"^^^ ■^«i/^^^"'*5)-^«2/i3K><3)c^^(?i?2 ^ Xi +X?) , (B.l) 



9192 



with Xi7 Xj being the two produced massive particles of mass rrii, rrij. Here fq^/A{xa^Q) 
is the distribution function of partons of type {qi = g, q, q), in the hadron of type A at a 
factorization scale Q. 

Taking the AB-c.m. system as the lab-system, the lab-momenta of the produced Xi 
and Xj are 

Pt = {Ei,PT, Pi cos 0i) , = {Ej, -pt, Pj cos Oj) , (B.2) 
where their transverse momenta are obviously just opposite 

Pt = Ptz = -PTj , (B.3) 
while their transverse energies Eti — VpfH-mf , Etj — -y/py -|- are used to define 



XTi = PTi = PT/ETi = -l / 1 - 



2 



/— ) r-'ii fi/^ii \l 2 ' 

V* V ^^Ti 



XTj = /3t,- = Pt/Etj = J 1 - • (B-4) 



Note that 



. ^ £;2 . ^ ^2 _ ^2 ^2 . ^ ^2 . ^ ^^^) (B_5) 



The rapidities and production angles of x^, Xj, in the lab-system, are related to their 
energies and momenta along the beam-axis of hadron A, (taken as the i-axis) by 

^22,, ^ Ei+ Pi COS Oj ^ ^^y. ^ E^+pj cose, 
Ei — Pi cos 9i ' Ej — Pj cos 9j 

The center-of-mass rapidity y of the x? Xj pair, and their respective rapidities y* in their 
own cm. frame, are defined as 

Vi^y + y* , yj^y + y*j , (b.t) 
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^y = y^-y^=yl-y* . (B.8) 

The fractional momenta of the the incoming partons are expressed in terms of their 
lab- momenta by (compare f|T| [3T1) ) 



gi = |(a;a,0,0,Xa) , g2 = |(a;f„ 0, 0, -Xfe) , q = q^ + q2 , 

g° = ^(^^a + Xb) = Ei + Ej , = ^(a^a - Xb) = {pi cos 9i + pj cos 9j) , (B.9) 
which lead to 

2 ^/s 

1 M - 

Xb = -[xne-y- + XTjC-y^] = -^e-y , (B.IO) 

s = = (gi + q2Y = XaXbS = + x'^j + 2xTiXTj cosh(Ay)] . (B.ll) 

Using this, s, Xa, Xb may be calculated in terms of the final particle rapidities yi, yj and 
their transverse momenta. From them, y is also obtained, and {y*, y*) through (lB.7p . 
The remaining Mandelstam invariants of the subprocesses satisfy 

i = {Pi - qiY = m\ - M{E* - p* cos 9*) = m\ - ^Mv^^e"^? = m\ - ^XaXne-y' 

= m] - M{E* - p* cos 9*) = m] - ^M^e^^ = m] - ^XbXrje^^ , (B. 12) 

u = {pj - q^f = m\ - M{E* + p* cos 9*) = ml - ^M^se^^ = m\ - ^XbXne^' 

= m]- M{E* + p* cos 9*)=m^^- ^Mv^e"^^* = - ^XaXTje-^' , (B. 13) 

J J 2 2 

s 

T = - = XaXb , (B.14) 

s 

where 9* describes Xi* production angle in the frame (the one being 71 — 9*). 

The energies of the two final particles in their c.m.-frame are 

s + mf — s + — mj 
E:= -^j^ , E*= -^7^, (B.15) 



their momentum is 



and their velocities 



P* = - m\ - m]f - Amlm^f , (B.16) 



Pi -P /Ei - — 



M2 + (mf - mh 



3' 

1 
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We also have 



cosr 



tanh y* 



tanh y* 
Pi 



sin e* 



p* 



Note that 



Xi 



Xj 



u 
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mf 



t — rrij 



u 



l + /3*cose* 
1 - P* cos 9* 
1 -/3;cosr 
1 + P* cos 9* 



p*cos9* 
Xj 



Pt 



u 



t Xi 



1 



ii + t Xi + 



1 ' 

+ M2 



E* 



cosh y* 



M2(l + x.)2 



Ti 



4(M2 + - m2)2;>^j 



X 



4(M2 + m- 



B.2 The basic distributions at LHC. 

Using fIB.ll) we define 

; On 



or 

da _M^ +mf- m] 

"^^ij 1 



dM'^dxtdy M2(l + Xi)2 
or 

da M 



dMHydpl 2s^p*^-pI ' 

where Sij describes the contribution to fIB.ip from all partons. Grouping tog 
gluon-gluon and qq parton contributions, we may write 

C — Q9 j_ cquark 
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with 

Sf, ^ 9{xa, Q)9{x,, Q)^^^^^^^ , (B.28) 



quark 



da{qq XiX^j) , . ^Ja{qq 



q{xa, Q)q{xb, Q) ,^ ' ^ + q{xa, Q)q{xb, Q) 



di di 



,(B.29) 



where q = u,d, s, c, b. 

For the numerical calculations presented here, we use as an example the MRST2003c 
code for quark and gluon structure functions |TT], taking the factorization scale as 



Q = Ell±^ . (B.30) 

As seen from above, the basic quantities needed s = M^, i,u. In case of 

(lR2l these are calculated from (iBTnj [RTTI [RT2l IRTI) . In case (ESSj), Eqs. (IRTqIIRTI 
[RT51 [RTHl [RTTI [RTHll must also be used. 



Starting from this basic distribution, always assuming > rrij, and imposing the 
cuts 

IVil < , \yj\ < Yj with Y, < Yj , (B.31) 

(in the numerical applications we take Yij = 2) we get the single variable distributions 
defined in the following subsections. 



B.2.1 The transverse energy and pr distributions 

^= jdmjdy,^S., , ^ = I dyj dy,^^S., , (B.32 
where s = is determined from (IB. 110 . xtj from (IB.SP and the integration limits are 



yjn^in = max \ In ( - — _ ) ; -Yj , , 



yjmax = min <! In ( - — ^^'^ ) ] Yj \ , (B.33) 



Vimax Uimir, 



2 - XTiC-y^ 

Xt 
XTj 



1 mf-m']\ (2-XTje~ 



min { Yi ■ cosh-i ( —(1 + ^ ; In ) |> . (B.34) 

.XTi S I \ XTi 



The XTi range in flB.32p would be 



2m,- mf - ^ ^ 

<XTi<l+ — , (B.35) 



'S s 

determined by the requirement that the middle constraint in flB.34p is meaningful. 
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B.2.2 The rapidity distribution 

Since the distribution has to be symmetric, we only consider the case of j/j > 0, for 
calculating the XTi-limits. Then in 

^ = J dxTi J dyj^^-^Sij , (B.36) 
the i/j-limits are given by (]B.33|) . while the limits for the xxi integration are 

XTimin •^Tmin,exp ^ "F^ ; 



9 9 

mf—m^ 



fi H — - — - ze ^ — v^i 1 

XTimax = Klin < — : — . ; — : 2e~^' > (B.37) 

I cosh?/i e2(^j±f») - 1 ' i ^ ^ 

with 

A, = 4[e2^^ + ^""^'"""^'^ 1 - e^^iJ^O)] , (B.38) 
provided ±yi > —Yj and (of course) yi > 0. 



B.2.3 The invariant mass distribution 

Using the cm. rapidity y defined above and fIB.lQp . one obtains 

da 



where the integration limits are 



(M2 + mf-m2) 
dxi / dy 7z 



X.ir. 



Xi' 



min \ Yi - ^In^i; Yj - ^In 



max 



r, - -In 

^ 2 



ln( 



mm 



1 + 13* M^{s + {m} - m^j)e 



M'^Xi — fTif + 



-2Yi 



1 - /?; ' M4e-2^^ + simf 



max 



2M2 

1 - (3* M^e-^^^ 



M2 



2K 



2M2 



M2 
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(B.39) 



(B.40) 



(B.41) 



with 



B.2.4 The angular distribution 

This is given by 



where the y integration hmits are as in flB.40p . while for the integration we have the 
hmits 

= max |Li; L2, L3; M!; MI^} , (B.44) 



with 



Ml = ix.K-m^ + se^^0±v/A3], (B.45) 



A3 = x'K' - ^? + se^^O' - 4sK' - mj)e'^% (B.46) 
1 

2x1 



Mf = ^[m^-mj + se'^^' i v^Ay, (B.47) 



A'3 = - m| + se^^O^ - 4s(m2 - m|)x,V^% (B.48) 

1 4rT7^ / 

= -[^-2K-m|) + v^], (B.51) 

X = (B.52) 



m7 mfimf - m^) 



4 

(1 - x'Y 1 - X 



A4 = 16 [ : ^\ \, (B.53) 
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Figure 1: Feynman diagrams for qq XiX?- 




s(TcV^) 



Figure 2: The s distributions in the SPSla model of [8]. 
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Figure 3: The s distributions in the SPS5 model of [8]. 
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Figure 4: The s distributions in the SPS6 model of [8]. 
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Figure 5: The s distributions in the SPSS model of [8]. 
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XTi 



Figure 9: The XTi distribution in SPSla of [8]. 




Figure 10: The xti distribution in SPS7 of [8]. 



35 




Figure 11: The Xi distribution in SPSla of [8]. 




Figure 12: The Xi distribution in SPS7 of [8]. 
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